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SUMMARY 


This  program  Is  being  conducted  for  the  purpose  of  develop- 
ing a predictive  model  of  heat  flow  and  solidification  for  alum- 
inum alloys  produced  under  the  high  cooling  rate  conditions 
achievable  In  rotary  atomization  and  laser  surface  melting.  It 
Is  a combined  experimental  and  theoretical  study  of  the  relation- 
ship between  the  Important  solidification  variables  (e.g.  cooling 
rate,  temperature  gradients.  Interface  shape  and  velocity,  super- 
cooling and  transformation  kinetics)  and  the  structure  and  micro- 
chemistry of  rapidly  solidified  aluminum  alloys. 

During  the  first  six  months  of  this  program  work\was  carried 
out  In  the  following  areas:  \ 

(a)  A computer  heat  flow  model  for  the  rapid  solidification 
of  metal  powders  was  developed. 

(b)  High  purity  A 1 - 4 . 5%  Cu  alloy  Ingots  were  prepared  and 
sent  to  United  Technology  Corporation  for  rotary  atomization.  Ini- 
tial analysis  of  the  structures  of  the  metal  powders  produced  were 
carried  out. 


(c)  Generalized  heat  flow  equations  and  computer  codes  were 
developed  for  one,  two  and  three  dimensional  heat  flow  during  rapid 
melting  and  subsequent  solidification  of  a surface  layer  subjected 
to  a high  intensity  heat  flux,  e.g.  that  generated  by  a continuous 
wave  C02  laser. 

(d)  Substrates  of  A 1 - 4 . 5%  Cu  alloy  were  specially  prepared 
and  sent  to  United  Technology  Corporation  for  surface  melting. 
Initial  characterization  of  the  rapidly  solidified  specimens  were 
carried  out  using  electron  microscope  techniques. 

Finally,  a 15KW  continuous  wave  C0?  laser  Is  presently  being 
installed  in  our  laboratories  and  should  become  available  for  use 
in  this  program  In  the  near  future. 
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I . ATOMIZATION  STUDIES 

These  studies  are  aimed  at  enhancing  our  current  limited 
understanding  of  rapid  solidification  phenomena  occuring  during 
atomization  of  alloy  powders  by  the  centrifugal  technique.  It 
Is  now  apparent  that  undercool i ng , nucleatlon  and  growth  pheno- 
mena during  rapid  solidification  can  not  always  be  explained  with 
available  theories  of  non-equilibrium  dendritic  solidification. 

For  example,  some  centrifugal ly  atomized  powders  show  microcrystal- 
line, non-dendr itic , structures  that  can  not  be  directly  correlated 
to  average  cooling  rate  during  solidification.  Other  phenomena, 
such  as  increased  solid  solubility,  repression  of  non-equilibrium 
second  phases  and  solute  enriched  dendritic  cores  are  influenced 
by  cooling  rate  prior  to,  as  well  as  during  solidification.  Full 
exploitation  of  these  observations  in  design  of  specific  high 
strength  alloy  compositions  for  rapid  solidification  processing 
requires  a much  better  understanding  of  the  relationship  between 
process  variables  - microstructure  and  microchemistry.  While  the 
specific  emphasis  of  the  work  in  this  program  is  on  aluminum  alloy 
powders  produced  by  centrifugal  (rotary)  atomization,  the  findings 
should  be  applicable  to  other  alloy  systems  and  atomization  processes 

In  the  first  six  months  of  this  program  a comprehensive  liter- 
ature survey  was  initiated  and  both  theoretical  and  experimental 
studies  were  carried  out  on  rapid  solidification  of  aluminum  alloy 
powders.  Initial  findings  from  this  work  are  presented  below. 

1 . Literature  Survey 

The  literature  survey  initiated  was  on  the  reported  effects 
of  rapid  solidification  on  the  microstructural  modi f i ca t i ons  of 
aluminum  alloys,  as  well  as  the  theoretical  models  available  that 
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relate  high  cooling  rates  to  nucleatlon  and  transformation  kin- 
etics during  crystalline  solidification.  Results  of  the  reported 
mlcrostructural  modifications  would  be  used  to  both  select  spe- 
cific aluminum  alloys  for  this  study  and  to  relate  and  compare 
the  effect  of  calculated  achievable  coollne  rates  In  rotary  atom- 
ization on  structure  and  microchemistry  with  other  rapid  solid- 
ification processes.  Reported  mlcrostructural  modifications  of 
aluminum  alloys  due  to  rapid  solidification  can  be  summarized  as 
follows: 

a)  Mlcrostructural  refinement,  manifested  as  smaller  grai.1. 
size  and  dendrite  arm  spaclnqs. 

b)  Extension  in  terminal  solid  solubility  of  the  primary 
a-Al  phase. 

c)  Morphological  changes  of  the  eutectic  or  the  primary 
phase. 

d)  Formation  of  non -equ i 1 ib r i urn  phases. 

e)  Coupled  eutectic  growth  at  off-eutectic  compositions. 

f)  Vacancy  super  saturation.  ' 

Reported  extensions  of  solid  solubility  and  non -equ  i 1 ibr  i unt 
second  phases  detected  in  binary  aluminum  alloys  are  summarized  in 
Tables  I and  II,  respectively.  A review  of  these  data  and  practical 
considerations  based  on  present  and  future  commercial  alloys  for 
rapid  solidification  applications  indicate  Al-Cu,  Al-Fe  and  Al-Si 
systems  to  be  most  promising  for  fundamental  atomization  studies. 

2.  Theoretical  Studies 


The  initial  studies  in  this  portion  of  the  program  were  aimed 
at  developing  a versitile  computer  heat  flow  model  for  solidifica- 
tion of  atomized  droplets  of  both  single  melting  point  metals  and 
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0.7 

923 

>6-9 

Ni 
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Formed  transiently  during  aging  following  a decomposition  kinetics 
similar  to  the  supersaturated  Al-Cu  solid  solutions. 
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alloys  that  freeze  over  a range  of  temperatures.  A generalized 
computer  code  was  developed  for  the  former  and  will  be  described 
herein,  while  work  on  the  alloy  system  Is  presently  underway  and 
will  be  reported  In  the  future. 

The  computer  heat  flow  model  developed  for  solidification  of 
pure  metals  and  alloys  that  freeze  at  a single  temperature  per- 
mits prediction  of  the  Important  solidification  variables  (e.g. 
cooling  rate,  Interface  velocity,  temperature  gradient,  etc.)  as 
a function  of  the  dimensionless  variables  governing  the  rate  of 
heat  extraction  from  metal  droplets.  A manuscript  based  on  this 
work  Is  in  preparation.  A short  summary  of  some  of  the  important 
findings  is  as  follows. 

In  general,  the  achievable  heat  transfer  coefficients  at  a 

c 

liquid  metal  d ropl et -env i ronment  Interface  are  limited  to  h < 10' 

2 (9) 

W/m  . K.  For  aluminum  metal  droplets  in  the  size  range  of  'v  lpm 
to  1000pm  this  translates  to  a limitation  on  the  range  of  Biot 
Numbers  of  interest,  10"^  > Bi  > 1.0. 

« 

hr 

M ■ T°  0)  . 

kL 

where  h Is  the  heat  transfer  coefficient  of  the  metal  droplet- 
environment  Interface,  r0  Is  the  radius  of  the  droplet  and  is 
the  conductivity  of  the  liquid  metal. 

Figure  1 shows  calculated  dimensionless  temperature  versus 
dimensionless  distance  in  a liquid  aluminum  metal  droplet  for  var- 
ious Biot  Numbers  at  the  instant  its  surface  reaches  the  melting 
point.  The  data  shows  that  for  Biot  Numbers  less  than  '\-0.01  there 
Is  no  significant  temperature  gradient  in  the  droplet  and  the  simple 

I 

] 
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Newtonian  cooling  expressions  can  be  used. 

An  Important  variable  effecting  supercooling  prior  to  nuc- 
leation  Is  the  cooling  rate  In  the  liquid  droplet.  A generalized 


expression  relating  the  Instantaneous  average  coollne  rate,  e 


a vg  * 


In  a liquid  metal  droplet  to  the  Blot  Number  and  dimensionless 
surface  temperature  was  derived: 


3 x B1  x 0 


SURFACE 


where : 


P - — 

ro 


= -3  { p1  <&>  d> 


T - T 


Tm  ' Tg 


t* 


T = temperature  in  the  droplet 

= temperature  of  the  environment 
Tm  = melting  point  of  the  droplet 

= thermal  diffusivity  of  the  liquid  metal 
t = time 

Computer  calculations  show  that  dimensionless  solidification 
interface  velocity  3p^/a(Fo)  increases  with  both  increasing  Biot 
Number  and  with  distance  solidified  from  the  droplet  surface  at 
a given  Biot  Number,  Figure  2. 

Examples  of  the  type  of  numbers  which  can  be  generated  from 
the  data  presented  in  equation  (2)  and  Figure  (2)  are  as  follows. 

A 50um  diameter  aluminum  droplet  convectively  cooled  with  a high 
velocity  gas  environment,  as  in  the  centrifugal  atomization  process 
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of  the  United  Technologies  Corporation,  may  achieve  a surface  con- 
ductance  (surface  heat  transfer  coefficient)  of  h ^5  x 1(j  W/m  .°K. 
The  calculated  average  cooling  rate  in  the  liquid  droplet  from 
equation  (2)  when  its  surface  reaches  the  melting  point  is  ^l.A  x 
10®  °K/sec.  The  solid-liquid  interface  velocities,  from  Figure  2, 

i 

at  radius  ratios  of  0.95  and  0.25  are  0.037  m/sec  and  0.32  m/sec, 
respectively. 

As  noted  above,  a new  heat  flow  model  is  presently  being  devel- 
oped for  solidification  of  alloys  that  freeze  over  a range  of  temp- 
eratures, hence  liquid  compositions. 

3.  Experimental  Studies 

Several  high  purity  ingots  of  Al-4.5%  Cu  alloy  weighing 
approximately  2 kilograms  were  prepared  in  our  laboratories  and 
sent  to  United  Technologies  Corporation.  Two  of  the  ingots  were 
remelted,  centri  fugally  atomized  and  solidified  by  forced  convec- 
tive cooling.  The  resulting  powders  are  being  subjected  to  micro- 
structural  analysis  using  various  electron  microscope  techniques. 
Figure  3 shows  the  three  distinct  types  of  structures  noted  to  date. 
The  micrograph  in  Figure  3(b)  shows  typical  dendritic  structures 
found  in  the  atomized  droplets.  Average  measured  dendrite  arm 
spacings  are  of  the  order  of  'v  2 pm.  The  micrograph  in  Figure  3(c) 
shows  a duplex  microstructure  - primary  spheroidal  solid  particles 
surrounded  by  a fine  dendritic  matrix.  This  structure  is  very 
similar  to  that  obtained  when  a metal  alloy  is  subjected  to  high 
rates  of  shear  during  initial  stages  of  solidification  - Rheocast 
structures  (10).  It  is  postulated  that  these  droplets  partially 
solidified  on  the  disk  of  the  rotary  atomizer  prior  to  flight  into 
the  forced  convective  cooling  environment.  On  the  other  hand,  the 
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powder  particle  In  Figure  3(a)  also  shows  a duplex  structure  which 
Is  similar  to  Figure  3(c)  yet.  It  is  significantly  different  In 
Its  response  to  etchants  and  its  general  appearance  In  the  Scan- 
ning Electron  Microscope. 

Detailed  structural  and  microchemical  analysis  of  these  struc- 
tures Is  presently  underway.  Techniques  are  being  developed  for 
the  preparation  of  electron  transparent  specimens  of  the  powder 
cross-sections  for  Scanning  and  Conventional  Transmission  Electron 
Microscope  Studies. 

Ingots  of  high  purity  Al-Fe  and  Al-Si  alloys  are  also  being 
prepared  in  our  laboratories  for  centrifugal  atomization  at  the 
United  Technologies  Corporation. 
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II . LASER  SURFACE  MELTING 

A corollary  theoretical  and  experimental  program  was  under- 
taken to  investigate  rapid  surface  melting  and  subsequent  solid- 
ification of  aluminum  alloy  substrates  subjected  to  high  intensity 
heat  fluxes  such  as  that  generated  by  continuous  wave  CO2  lasers. 

The  two  features  of  rapid  surface  melting  which  make  it  an  attrac- 
tive tool  to  correlate  rapidly  solidified  structures  to  solidifi- 
cation variables  are: 

(a)  Theoretical  heat  flow  calculations  can  be  carried  out 
without  recourse  to  arbitrary  assumption  of  a heat  transfer 
coefficient  between  the  liquid  and  the  substrate  below  - 
there  is  intimate  contact  between  the  two. 

(b)  It  separates  growth  phenomena  from  nucleation  events 
in  rapid  solidification  processing  - special  substrates  can 
be  prepared  for  epitaxial  growth  studies. 

Our  investigations  are  aimed  at  coupling  computed  heat  flow 
conditions  with  experimentally  observed  microstructures  in  the 

1 

aluminum  alloys  used  for  the  centrifugal  atomization  studies.  Here- 
tofore, laser  surface  melting  of  our  aluminum  substrates  was  car- 
ried out  at  the  United  Technologies  Corporation.  Recently,  a 15KW 
continuous  wave  CO^  laser  was  acquired  by  the  University  of  Illin- 
ois and  is  presently  being  install ed  in  our  laboratories.  It  is 
anticipated  that  specimens  would  also  be  surface  melted  with  this 
new  in-house  facility. 

1 . Heat  FI ow  Cal cu  1 a t i ons 

Generalized  heat  flow  equations  and  computer  codes  were  de- 
veloped for  prediction  of  the  important  rapid  melting  and  subsequent 
solidification  variables  of  surface  layers  subjected  to  high  intensity 


heat  fluxes.  Specific  work  carried  out  In  the  first  six  months 
of  this  program  is  summarized  below. 
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( a ) One-Dimensional  Hea t Flow 

Numerical  calculations  of  a one-dimensional  heat  flow  model 
developed  earlier  were  completed  and  published  during  this  contract 

I 

period.  The  original  manuscript  submitted  for  publication  prior 
to  this  contract  was  revised  to  include  these  recent  calculations. 
The  published  manuscript  is  attached  hereto  as  an  Appendix. 

Shortcomings  of  this  work  included  assumptions  of  uniform 
absorbed  heat  flux  and  one-dimensional  heat  flow  - it  is  assumed 
that  the  laser  spot  diameter  is  large  in  comparison  with  the  heat 
affected  zone  in  the  substrate. 

( b ) Two-Dimensi onal  He at  Flow 

In  this  study  we  have  developed  generalized  expressions  for 
the  finite  difference  coefficients  of  the  energy  conservation 
equation  in  a curvilinear  coordinate  system.  We  have  considered 
the  rapid  melting  and  subsequent  solidification  of  the  surface  layer 
of  a semi -inf ini te  aluminum  substra te , . i ni tial ly  at  room  tempera- 
ture, subjected  to  a high  intensity  heat  flux  over  a circular 
region  on  its  bonding  surface.  Both  uniform  and  Gaussian  distri- 
butions of  the  absorbed  heat  flux  were  considered.  However,  the 
equations  and  computer  codes  developed  can  treat  any  arbitrary 
heat  flux  over  the  circular  region. 

Computer  calculations  are  presently  underway  for  a range  of 
absorbed  heat  flux  and  the  results  of  the  computations  are  being 
compared  with  the  one-dimensional  heat  flow  model.  Figure  4 shows 
the  calculated  fractional  melt  depths  (liquid-solid  interface  pos- 
itions) at  various  times  during  surface  melting  of  an  aluminum 
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substrate  subjected  to  a uniform  absorbed  heat  flux  of  5 x 10 
W.m  . The  two  sets  of  curves  are  for  different  laser  spot  radii, 
190pm  and  380pm.  Figure  5 shows  similar  data  obtained  for  two 
different  absorbed  heat  fluxes.  The  data  in  Figures  4 and  5 In- 
dicate that  some  of  the  general  relationships  developed  between 
absorbed  heat  flux,  melt  depth  and  time  in  the  one-dimensional 
model  (see  Appendix)  may  also  hold  for  the  two-dimensional  trans- 
ient heat  flow  model. 

It  is  anticipated  that  data  generated  in  the  next  two  months 
would  permit  the  development  of  specific  relationships  between  ab- 
sorbed heat  flux  and  isotherm  shapes,  isotherm  velocities  and  temp- 
erature gradients.  Experiments  are  also  planned  in  which  the  con- 
ditions for  these  two-dimensional  transient  heat  flow  calculations 
can  be  simulated  by  laser  radiation  of  stationary  aluminum  substrates. 

( c ) Three-Dimensional  Heat  Flow 

In  this  work  we  have  addressed  the  problem  of  rapid  surface 
melting  and  subsequent  solidification  of  a semi-infinite  solid, 
moving  with  constant  vel oc i ty,  s ub jected  to  a high  intensity  heat 
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flux  over  a circular  region  of  its  bonding  surface.  Temperature 
profiles  in  the  molten  region  and  the  adjacent  heat  affected  zone, 
as  well  as  the  important  melting  and  solidification  variables  are 
to  be  determined. 

Generalized  expressions  for  the  coefficients  of  the  finite 
difference  equations  governing  heat  transfer  in  discretized  spatial 
domains  in  moving  curvilinear  systems  have  been  developed  along 
with  a computer  code  for  numerical  solutions.  Initial  data  ob- 
tained  show  that  the  method  used  is  accurate  and  does  lend  itself 
to  solution  of  this  multi-dimensional  phase  change  problem.  It 
is  anticipated  that  this  Investigation  would  be  completed  in  the 
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next  six  months  of  this  program. 

Solution  of  the  problem  described  above  should  accurately 
simulate  the  thermal  field  geometries  encountered  in  laser  surface 
melting  of  moving  substrates.  This  Is  the  technique  presently 
used  in  our  experimental  work.  Finally,  the  general  expressions 
and  solution  techniques  developed  should  also  be  applicable  to 
calculation  of  thermal  fields  and  solidification  parameters  In 
the  metal  pool  of  certain  types  of  welding  processes. 

2.  Experimental  Stud ies 

Unidirectional  ingots  of  Al-4.5%  Cu  alloy  were  cast  and  sec- 
tioned into  appropriate  sizes  for  laser  surface  melting.  These 
were  surface  melted  at  United  Technologies  Corporation  under  con- 
ditions predicted  from  the  heat  flow  studies  to  yield  a range  of 
melt  depths  and  solidification  conditions.  Preliminary  examina- 
tion of  the  structures  obtained  via  SEM  and  Conventional  TEM  tech- 
niques has  been  completed.  Figure  6 shows  representative  TEM 
photomicrographs  of  one  laser  melted  specimen.  As  anticipated,  the 
structures  are  very  fine  and  are  aligned  in  the  heat  flow  direction. 
Work  is  now  underway  to  determine  composition  profiles  across  the 
fine  segregate  spacings  via  STEM.  Correlation  of  structures  and 
microchemistries  to  calculated  heat  flow  conditions  and  solidifi- 
cation - segregation  models  would  be  an  important  aspect  of  this 
work  in  the  future. 

Presently,  special  single  crystal  specimens  of  aluminum  alloys 


are  being  grown  in  our  laboratories  for  laser  surface  melting  stud- 
ies. It  is  expected  that  epitaxial  growth  studies  could  be  better 
carried  out  on  these  substrates. 
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Dimensionless  Temperature 


Figure  1 Calculated  dimensionless  temperature  versus  dimension- 
less distance  in  liquid  droplets  of  aluminum  for  dif- 
ferent Biot  Numbers. 


Biol  Number  (hr0/kL) 

Figure  2,  Dimensionless  solid-liquid  interface  velocity  versus 
Biot  Number  during  solidification  of  aluminum  powders 
at  different  interface  positions. 


Figure  3 Scanning  Electron  Micrographs  of  three  distinct 

microstructures  observed  in  Al-4.5*  Cu  alloy  pow- 
ders produced  by  the  centrifugal  atomization  tech- 
nique. 
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APPENDIX 


Rapid  Melting  and  Solidification  of  a Surface 
. Layer 

S.  C.  HSU.  S.  CHAKRA VORTY.  AND  R.  MEHRABLAN 

A one-dimensional  computer  heat  flow  model  is  used  to  investigate  the  effect  of  high  in- 
tensity heat  fluxes,  e.g.  those  achieved  via  continuous  CO-  laser  radiation,  on  the  im- 
portant surface  layer  melting  and  subsequent  solidification  variables  of  three  substrate 
materials:  aluminum,  iron,  and  nickel.  Temperature  profiles  rs  time,  melting,  and  soli- 
dification interface  velocities,  heating,  and  cooling  rates  in  the  surface  layers  of  the 
three  metals  are  calculated.  Results  are  presented  in  a general  form  to  permit  deter- 
mination of  these  variables  for  large  ranges  of  absorbed  heat  fluxes  and  times.  General 
trends  established  show  that  temperature  gradients  in  the  liquid  and  solid  phases  and  in- 
terface velocities  are  directly  proportional  to  the  absorbed  heat  flux,  whereas  melt  depth 
is  inversely  proportional  to  the  absorbed  heat  flux.  Average  cooling  rates  comparable  to 
splat  cooling  can  be  achieved  by  increasing  the  heat  flux  and  reducing  the  dwell  time  of 
the  incident  radiation.  An  order  of  magnitude  increase  in  the  absorbed  heat  flux  results 
in  a corresponding  two  orders  of  magnitude  increase  in  average  cooling  rates  in  the 
liquid  during  solidification  of  crystalline  and  noncrystalline  structures. 


I INTRODUCTION 

The  effect  of  high  cooling  rates.  10'  to  10‘"  K s on 
the  structure  of  solid  materials  formed  from  their 
melts  has  received  increased  attention  ever  since 
Duwez  and  his  coworkers  published  results  of  their 
investigations  on  splat  cooling.1  ,s  A large  number  of 
innovative  batch  and  continuous  techniques  for  pro- 
duction of  laboratory  and  commercial  quantities  of 
rapidly  solidified  material  have  since  been  developed. 

In  general,  these  techniques  involve  quenching  of  a 
liquid  on  a high  conductivity  substrate.  Determina- 
tion of  exact  cooling  rates  during  solidification  of  crys- 
talline and  noncrystalline  structures  in  these  pro- 
cesses has  required  estimates  of  heat  transfer  coeffi- 
cients between  the  liquid  and  the  substrate.  Two  note- 
worthy attempts  have  been  made  to  directly  measure 
heat  transfer  coefficients  during  splat  cooling. How- 
ever. it  is  questionable  whether  these  values  can  be 
used  in  calculation  of  cooling  rates  in  other  rapid  soli- 
dification techniques  or  even  in  splat  cooling  in  general 
because  of  different  processing  conditions,  alloy  com- 
positions and  substrate  materials. 

The  recent  availability  of  high  power  continuous 
wave.  CO-  laser  has  led  to  the  development  of  a rapid 
surface  meltingand  solidification  technique'  in  which 
the  bulk  (semiinfinite)  substrate  in  intimate  contact 
with  the  molten  layer  acts  as  the  quenching  medium - 
an  infinite  heat  transfer  coefficient  can  be  assumed. 
This  paper  describes  a theoretical  study  carried  out 
to  investigate  the  effect  of  high  intensity  radiation  on 
the  important  surface  layer  melting  and  subsequent 
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solidification  variables  of  three  substrate  materials: 
aluminum,  iron,  and  nickel. 

Several  studies  have  previously  been  carried  out  to 
quantitatively  describe  heat  flow  in  a semiinfinite 
substrate  material  subjected  to  high  intensity  radia- 
tion. An  analog  computer  program  was  developed  by 
Cohen1'  to  investigate  the  temperature  distributions  in 
the  molten  surface  layer  and  the  solid  substrate  dur- 
ing melting  of  a semiinfinite  solid  subjected  to  a step- 
function  heat  input.  The  thermal  profiles  thus  gen- 
erated were  subsequently  used  by  other  investigators 
to  estimate  average  cooling  rates  during  solidification 
of  noncrystalline  structures.5  Other  researchers7 
have  used  the  point  source  heat  flow  equations  pro- 
posed by  Rosenthal"  for  bead-on  plate  welding.  One 
difficulty  with  all  theoretical  calculations,  including 
that  presented  here,  is  that  not  all  incident  radiation 
is  absorbed  by  the  substrate- refle  ivity  varies  with 
power  density,  substrate  material,  surface  condition 
(e.g..  oxide  layer),  and  is  different  for  primary  phases 
and  inclusions.  Nevertheless,  it  was  believed  that  a 
detailed  study  of  the  thermal  characteristics  of  pure 
substrate  materials  would  yield  useful  relationships 
between  the  absorbed  heat  flux  and  the  important  melt- 
ing and  solidification  variables  influencing  structure. 

In  practice,  dwell  time  of  an  incident  radiation  can 
be  calculated  from  the  scanning  speed  and  spot  diame- 
ter of  the  incident  beam.5,7  The  one-dimensional  heat 
flow  model  used  in  this  study  assumes  that  the  spot 
diameter  is  large  in  comparison  with  the  melt  depth. 
Furthermore,  it  is  assumed  that  there  is  no  convec- 
tion in  the  molten  surface  layer.  The  effect  of  con- 
vection on  melt  depth  was  determined  by  arbitrarily- 
increasing  the  conductivity  of  the  liquid.  In  the  ex- 
treme. the  conductivity  of  a well  mixed  liquid  tends  to 
infinity  and  analytical  equations  developed  by  Landau" 
for  complete  removal  of  the  melt  could  be  applied.  In 
lieu  of  exact  analytical  solutions  for  the  problem  at 
hand,  approximate  analytical  techniques  such  as  the 
•'approximate  integral  technique"  developed  by  Good- 
man10 can  be  used.  Melt  depth  calculations  based  on 
this  latter  approach  were  also  carried  out  in  this 
study  and  are  compared  to  those  calculated  via  the 
numerical  integration  technique. 
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TabU  I.  Summary  of  Properties  and  Numerically  Calculated  Constants  ' _ 


Meliing 

Temperature.  7'^/ 

K 

Boiling 

Temperature,  Tv 

K 

Element 

Solid  Density. 

Kg  til1 

l iquid  Density,  *>/ 
Kg  m 1 

Conductivity 
of Solid.  ks 

J s' 1 m 1 1C  1 

Conductivity 
of  Liquid,  A/ 

J s'1  in  1 K” 

Aluminum 

2700 

2390 

933 

272.1 

228* 

108+ 

Iron 

7860 

7030 

1808 

3273 

42* 

44+ 

Nickel 

8900 

7900 

1728 

3005 

74* 

43 

Specific  Heat 

Specific  Heat 

Thermal  Dillusivity 

Themial  Dilfusivity 

Latent  Heat 

1.  + S 

'•  max 

Element 

of  Solid.  (' 

I Kg'  C' 

of  l iquid.  ( 

J Kg  ' K4 

of  Solid,  u, 

of  l iquid.  a/ 

of  fusion.  // 

m*  s'1 

m>  »-■ 

JKg1 

Aluminum 

1048* 

1086* 

8.1  X Iff1 

4 2 X Iff* 

3.95  X IO’5 

9.3 

Iron 

691* 

628t 

7 7 X I0‘ 

IOX  Iff5 

2.72  X 10s 

122 

Nickel 

556* 

65bt 

1.5  X Iff1 

8.3  X 10* 

2.99  X I05 

21.6 

W't,,, 

<b/'max 

P/M 

<rF-  rAf>c/)/ 

v/S/2v^7vs(rM-r(,) 

J m 1 s 1 1 

J m 1 s 1 * 

J 111 J s 1 1 

// 

J in  - s 1 ’ 

Aluminum 

i ox  io’ 

4.21  x io’ 

5 tux  10’ 

0 53 

4.92 

1.43  X 10’ 

Iron 

202  x io1 

4.07  X 10’ 

4.26  X 10’ 

0 23 

3.38 

202  X 10’ 

Nickel 

:.4.ix  io7 

4.27  X 10’ 

4.41  X 10’ 

0.33 

2.80 

2.43  X 10’ 

•Aveiaged  from  298  K to  M P 
t Averaged  from  M P ro  B P 


n.  HEATING  OF  THE  SOLID  III.  APPROXIMATE  INTEGRAL  TECHNIQUE 

Integration  of  the  one-dimensional  transient  heat 
conduction  equation  for  the  case  of  a semiinfinite  solid 
at  its  melting  point  subjected  to  a constant  heat  flux 
at  its  surface  leads  to  a heat  balance  integral  given 
by  Goodman:10  Solution  of  this  integral  with  the  ap- 
propriate boundary  conditions  at  the  surface  of  the 
metal  and  the  liquid-solid  interface  is  facilitated  by  an 
a priori  assumption  of  a temperature  distribution  in 
the  liquid.  In  this  work  the  temperature  distribution 
in  the  liquid  was  represented  by  a quadratic  equation 
and  the  analytical  solutions  to  the  heat  balance  integral 
derived  by  Evans  cl  al'z  for  this  case  were  used  to 
generate  curves  of  maximum  melt  depth.  L max . rs  ab- 
sorbed heat  flux.  The  maximum  melt  depth  throughout 
this  paper  refers  to  the  melt  depth  achieved  for  a 
given  absorbed  heat  flux  when  the  surface  of  the  ma- 
terial reaches  its  vaporization  temperature.  Tv,  and 
the  heat  flux  is  removed.* 

•See  Nomenclature  for  definition  of  all  notations. 

Figure  1 shows  calculated  maximum  melt  depths  VS 
absorbed  heat  flux  for  the  three  elements.  The  more 
precise  calculations  using  the  computer  model  dis- 
cussed in  the  next  section  are  also  shown  in  the  same 
figure.  The  procedure  outlined  above,  while  providing 
useful  order  of  magnitude  estimates  of  melt  depths, 
has  severe  shortcomings.  It  does  not  permit  deter- 
mination of  temperature  profiles  in  the  solid;  there- 
fore. it  cannot  be  used  for  calculation  of  thermal  pro- 
files during  solidification  of  the  molten  surface  layer 
once  the  heat  flux  is  removed. 

IV.  COMPUTER  SOLUTION 

The  numerical  integration,  finite  difference  method, 
developed  by  Murray  and  Landis10  was  used  for  both 
the  melting  and  subsequent  solidification  of  a surface 
layer  subjected  to  a constant  heat  flux  for  a given 
period  of  time.  Initial  conditions  were  established  us- 
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The  temperature  distribution  in  a semiinfinite  solid 
subjected  to  a constant  heat  flux  at  its  surface  is 
readily  determined  from  equations  given  in  Carslaw 
and  .laeger.11  These  equations  were  used  to  establish 
the  initial  temperature  distributions  rs  time  in  the 
different  substrates  prior  to  surface  melting.  The 
product  of  where  <i  and  /,„  are  absorbed  heat 

flux  and  time  for  the  surface  of  the  solid  to  reach  its 
melting  point,  respectively,  is  a constant  for  a given 
material.  The  calculated  values  of  this  constant  and 
all  other  properties  of  the  three  elements  (aluminum, 
nickel  and  iron)  used  in  this  study  are  listed  in  Table  I. 

Data  thus  generated  were  used  both  in  the  deter- 
mination of  melt  depths  rs  time  via  the  approximate 
integral  technique  and  to  establish  initial  temperature 
profiles  for  the  numerical  calculations. 


ABSORBED  MCA r FLU*.  U m” 


Fig.  1— Maximum  melt  depth  rs  absorbed  heat  flux  for  alu- 
minum, iron  and  nickel  calculated  using  the  approximate  in- 
tegral and  numerical  techniques. 


mu  iltta  generated  from  analytical  expressions  dis- 
cussed earlier  (or  the  single  phase  solid  material 
Natural  convection  due  to  buoyancy  forces  in  the  liquid 
, ami  fluid  motion  due  to  density  differences  in  the  two 
phases  were  ignored  * Pure  metal  substrates 

•|hc  effect  V»l  convection  »*n  melt  vkpth  ykteinuned  in  * ic*  uint*iir  vaUuIj 
t»oo»  I he  te>ult\  v'l  ihew  v*kuiattom  will  he  t|UAitu<(t'rU  Jen.  nbtd 

(aluminum.  iron,  and  nickel!  were  considered.  A de- 
scription of  the  Murray  and  l.andis  variable -space 
network  technique  and  the  boundary  and  initial  con- 
ditions used  is  presented  below,  followed  by  the 
numerical  results  obtained  using  CDC  CYBER  175 
digital  computer  Thermal  and  physical  properties 
were  assumed  constant  but  different  for  the  liquid 
and  solid  phases.  Data  used  are  listed  in  Table  1. 


1'  Differential  and  Difference  Equations  for  a 
Moving  Boundary 

The  differential  equations  for  one  dimensional  heat 
flow  with  constant  but  different  thermal  properties  m 
the  liquid  and  solid  regions  are  coupled  at  the  liquid- 
solid  interface  by 

*-(•*)  (>HR  1 

where  K is  the  velocity  of  the  liquid-solid  interface; 

It  is  the  heat  of  fusion  for  the  pure  metal  (negative 
for  melting  and  positive  for  solidification);  and  p is 
the  density  of  the  liquid  for  melting  and  density  of  the 
solid  for  solidification. 

The  liquid  ami  solid  pluses  were  divided  into  equally 
sized  space  increments  which  changed  in  size  as  the 
liquid-solid  interface  moved  forward  during  melting 
and  reversed  direction  during  solidification  of  the 
molten  layer  sometime  after  the  heat  (lux  was  re- 
moved from  the  surface.  The  three  point  finite  dif- 
ference approximations  of  the  one  dimensional  heat 
flow  equations  and  Eq  1 ( derived  by  Murray  and 
l.andis"  were  used 

In  all  the  calculations  the  distance  from  the  sub- 
strate surface  to  the  last  nodal  point  in  the  solid  was 
large  enough  so  that  the  temperature  of  this  nodal 
(xunt  remained  at  the  initial  temperature  of  the  solid. 

This  distance  was  determined  by  comparing  the 
calculated  temperature  of  the  last  node  to  l'„.  Finally, 
trial  ami  error  showed  that  an  efficient  utilization  of 
computer  time  without  loss  of  accuracy  required  10 
increments  in  the  liquid  and  20  increments  in  the 
solid  with  the  liquid- solid  interface  located  at  node 
No.  11. 


2)  Boundary  and  Initial  Conditions 

The  boundary  conditions  at  the  surface  and  last 
node  were  determined  by  carrying  out  an  energy 
balance  across  half  space  increments  at  these  loca- 
tions. 
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Equations  '2|  amt  (3|  are  for  the  surface  and  the  last 
node,  respectively . V is  the  total  number  of  males,  < 
is  the  number  of  space  increments  m the  liquid  re 
gton.  I is  the  location  of  the  liquid-solid  interface 
at  r * 1.  t>  I'  I',,,.  and  E is  the  location  of  the  last 
node  in  the  solid. 

The  additional  boundary  condition  at  the  liquid  solid 
interface.  V t . is 

/m -o  HI 

Heat  loss  due  to  radiation  and  convection  at  the 
metal  air  interface  was  included  in  the  calculations  by 
assuming  an  effective  heat  transfer  coefficient , h at 
this  interface  an  order  of  magnitude  larger  than  that 
calculated  from  radiation  only  The  term 

h(/’  /’„)  was  incorporated  m Eq  2 above.  The 

magnitude  of  this  heat  loss  is  negligible  compared  to 
the  absorbed  heat  fluxes  used  in  our  calculations 
Hence,  the  heat  loss  due  to  convection  and  radiation 
had  a negligible  effect  on  the  calculated  temperature 
distributions 

When  only  one  phase  is  present  prior  to  the  melting 
of  the  solid  or  at  the  end  of  solidification  \t  O', 
terms  in  the  finite  difference  equations  of  Murray 
and  l.andis"  and  Eq  2|  above  become  undefined 
Therefore,  the  melting  problem  must  be  started  with 
a small  initial  value  of  l.  and  a starting  temperature 
profile  in  this  liquid  region  must  be  specified  Simi- 
larly. the  solution  during  the  solidification  half  of  the 
program  must  be  stopped  before  l 0. 

A starting  procedure  similar  to  that  suggested  by 
lleitz  and  Westwater**  was  especially  developed  for 
the  specific  problem  of  this  investigation.  As  previ- 
ously noted  temperature  distributions  for  the  initial 
heating  of  the  solid  were  determined  from  equations 
given  by  Carslaw  and  Jaeger."  The  solution  to  these 
equations  was  continued  until  the  surface  of  the  solid 
reached  a temperature  above  the  melting  temperature 
of  the  metal  (the  solid  was  superheated'  Then  tin 
sensible  heat  of  the  superheated  solid  region.  \ V, , 
was  instantly  converted  to  heat  of  fusion,  resulting  in 
an  initial  melt  depth.  Y Y Y,  of  uniform  tempera 
ture.  /'  T,\t  . The  mathematical  equivalent  of  this 
statement  is 

Ps-l>s  f X\l'  15  i 

Equation  [5|  was  directly  integrated  after  substitut- 
ing expressions  given  in  Ref  It  The  resulting  equa 
tion  is. 

PjHY..  p«C©,ir„  /'u  Y i ♦ q/  ( I 4i  cue  1 ^ 

V 2s  V 

b 

Tabulated  values  of 

a-'  Y> 

4i  cue 

2 s 

are  given  in  Ref  15. 

Figure  2 shows  an  example  of  initial  temperature 
distributions  established  in  tins  way  m an  aluminum 
substrate  subjected  to  an  absorbed  heat  flux  of  5 v It' 

W m"  . Curves  2 and  3 m this  figure  are  for  the  super 
heated  solid  and  after  conversion  to  an  initial  melt 
depth,  respectively 
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Kin  2— Temperature  distribution  in  an  aluminum  substrate  at 
different  times  during  heating  of  the  initial  sol ul , melting  and 
subsequent  solidification  for  an  absorbed  surface  heat  flux  of 
T x 10’  \V, m2 . Temperature  distributions  at  successive  times 
are  denoted  by  numbers  1 through  8.  Last  nodal  point  in  the 
solid  is  at  a distance  of  ~600  pm. 


The  value  of  initial  melt  depth.  X..  used  to  start  the 
differential  equations  was  consistently  small  (between 
0.03  and  0.1  of  the  melt  depth).  These  starting  thick- 
nesses were  selected  to  both  preserve  accuracy  of  the 
temperature  profiles  generated  as  well  as  keep  solu- 
tion times  required  reasonable.  Computations  were 
carried  out  with  a time  interval  continuously  determined 
front  <r  A/ 1'  L~  - 1 4. 

The  method  used  to  generate  initial  melt  depths  of 
constant  temperature  T = T \/  results  in  an  initial 
negative  interface  velocity.  However,  the  melt  velo- 
cities rapidly  change  sign  and  increase  to  a maxi- 
mum value  over  a length  of  time  of  approximately 
l O'"  s. 


3)  Sequence  of  Numerical  Calculations 

The  sequence  of  calculations  carried  out  were  as 
follows:  Temperature  profiles  rs  time  were  estab- 
lished during  heating  of  the  solid  until  its  surface 
reached  and  exceeded  the  melting  point  of  the  material. 
The  superheated  solid  region.  X,.  was  converted  to 
initial  melt  depth.  X3.  using  Eq.  |5|.  Melting  was  con- 
tinued until  a desired  melt  depth  was  reached  or  until 
the  surface  of  the  liquid  layer  reached  the  vaporiza- 
tion temperature  of  the  metal  (see  curve  No.  5 in  Fig. 
2).  The  heat  flux  was  removed,  q - 0 in  Eq.  [2].  and 
calculations  were  continued  until  the  interface  e lo- 
city.  R.  became  negligible*  heat  flux  in  the  liquid  and 

•The  interface  velocity  decreased  to  a value  of  approximatelv  thtec  orders  ot 
injjamude  lew  than  the  average  for  a pven  power  input 

solid  at  the  liquid-solid  interface  became  equal.  The 
melt  depth.  L (f.max  when  the  heat  flux  was  re  moved 
at  Ts  - Tv)  was  thus  established  (see  curve  No.  6 in 
Fig.  2).  The  sign  of  the  heat  of  fusion.  H.  in  Eq.  l| 
was  changed  vj  commence  solidification  tsee  curves 


No.  7 and  8 in  Fig.  2).  Distance  solidified  is  denoted 
as  t it  is  the  distance  from  the  melt  depth,  /. . to  the 
location  of  the  liquid-solid  interface.  Time  in  all  the 
data  presented  herein  is  total  time  elapsed  from  initial 
application  of  the  heat  flux  to  the  solid  surface. 


V.  RESULTS  OF  NUMERICAL  CALCULATIONS 

The  procedure  described  in  the  previous  section  was 
used  to  obtain  temperature  profiles  rs  time,  melting, 
and  solidification  interface  velocities  and  heating  and 
cooling  rates  in  the  surface  layers  of  the  three  metals 
(aluminum,  iron,  and  nickel).  Results  of  this  work  are 
presented  in  a general  enough  form  to  permit  deter- 
mination of  most  of  these  variables  for  large  ranges 
of  absorbed  heat  fluxes  and  times. 

Calculated  melt  depths  rs  total  time  for  several  ab- 
sorbed heat  fluxes  are  shown  in  Fig.  3.  The  arrows 
in  the  figure  denote  times  at  which  the  surface  of  each 
material  reaches  its  vaporization  temperature  and  the 
heat  flux.  </.  is  removed.  As  noted  earlier,  melting 
continues  a while  longer  until  a maximum  melt  depth. 

I- mux  - 's  reached.  Maximum  melt  depths  rs  absorbed 
heat  flux  calculated  via  numerical  techniques  are 
plotted  in  Fig.  1 for  comparison  with  those  predicted 
using  the  integral  profile  technique.  Melt  depth  during 
the  melting  and  subsequent  solidification  periods  can 
also  be  expressed  in  dimensionless  form.  psHL 
Figure  4 shows  this  function  plotted  rs  dimensionless 
time,  t l,„.  for  the  three  elements.  The  maxima  in  the 
bell-shaped  curves  correspond  to  maximum  melt 
depths.  Lmax. 

The  ratios  of  melt  depth  to  maximum  melt  depth. 

L L max.  are  plotted  rs  </v  / for  the  three  materials  in 
Fig.  5.  On  the  right-hand  side  of  Fig.  5.  data  are  pre- 
sented to  permit  determination  of  distance  from  the 
interface  to  the  last  nodal  point  in  the  solid.  S.  at  any 
time  during  melting  or  solidification.  The  total  dis- 
tance from  the  surface  of  the  metal  to  the  last  nodal 
point  in  the  solid  is  E = S + L:  it  is  a constant  at  a 
given  absorbed  heat  flux.  F or  example,  an  aluminum 
substrate  absorbing  a heat  flux  of  i/  5 x 10'’  w nr' 

10" 
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Fig.  Melt  depth  » total  time  for  different  absorbed  heat 
flu'cs  obtained  via  'he  numerical  technique.  Arrows  indicate 
times  at  which  a surface  reaches  its  vaporization  tempera- 
ture and  heat  flux  is  removed. 
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Kig.  4— Dimensionless  melt  depth  psHL/qtm  during  melting 
and  solidification  vs  dimensionless  time, 
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Kig.  5— Dimensionless  melt  depth,  L/Lmax,  and  dimension- 
less distance  from  the  liquid-solid  Interface  to  the  last  nodal 
point  in  the  solid  vs  </V  I for  aluminum,  iron  and  nickel  sub- 
strates. 


will  achieve  a maximum  melt  depth  of  approximately 
65  pm  at  a time  of  / max  “ 1.08  x 10"1  s.  The  value  of 
S at  tnis  time,  t/Lmax  - 1.0,  is  obtained  from  Fig.  5; 
Syntax  = 8-3,  S 535  pm.  Therefore,  the  distance 
from  the  surface  of  the  metal  to  the  last  nodal  point 
is  E ^ 600  pm.  Using  the  curves  and  the  right-hand 
data  presented  in  Fig.  5,  the  melt  depth  and  the  heat 
affected  zone  in  the  solid,  S,  can  be  determined  any 
time  during  the  melting  and  the  solidification  period. 

The  information  generated  above  permits  deter- 
mination of  the  location  of  each  nodal  point  during 
melting  and  solidification.  The  temperature  profiles 
in  the  two  phases  throughout  the  melting  process  and 
during  solidification  after  reaching  a maximum  melt 
depth  can  thus  be  established  from  Figs.  6,  7,  and  8. 
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As  expected,  there  is  significant  loss  of  superheat 
from  the  liquid  region  during  the  time  period  /„  £ / 

- / mux.  On  fl'e  other  hand,  the  solid  portion  of  the 
substrate  heats  up  continuously  during  the  melting  and 
solidification  cycles. 

Interface  velocity,  li,  is  directly  proportional  to  the 
absorbed  heat  flux.  R/ q vs  fractional  distance  melted 
or  solidified,  X/ L max,  for  the  three  substrates  are 
shown  in  Fig.  9.  X denotes  melt  depth,  L,  during 
melting  and  distance  solidified,  e,  during  solidification. 
The  change  in  temperature  with  time  in  the  liquid  at 


Kig.  8— Dimensionless  temperature  at  each  nodal  point  in  Ihe 
liquid  and  the  solid  regions  as  a function  of  dimensionless 
time,  ///<«,  during  melting  and  solidification  of  the  surface 
layer  of  an  aluminum  substrate.  1 = 1 is  the  nodal  point 
located  at  the  surface  and  t = 31  is  the  last  nodal  point  In  the 
solid. 


Kig.  7— Dimensionless  temperatures  at  each  nodal  point  in 
the  liquid  regions  as  a function  of  dimensionless  time,  t/tm, 
during  melting  and  solidification  of  surface  layers  of  iron  and 
nickel  substrates.  » l is  Ihe  nodal  point  located  at  the  sur- 
face. 
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Ktg.  s— Dimensionless  lemperuture  at  several  nodal  points 
lit  ll\i>  solid  region  iih  ii  function  i>l  dimensionless  Mini',  1,1)11. 
during  melting  mill  solidification  Hi  surface  layers  of  linn 
mill  nickel  substrates.  i 1 1 In  tin’  liquid-solid  InU'i'lin't', 

/'  l',l| , i III  In  lilt'  liinl  nodal  point  In  l h«<  solid. 


Kill.  :)  I In'  I'tilln  nl  Interlace  velocity  In  absorbed  liniil  (lux, 
K ,/,  tin  n luiit-l Imi  nl  fractional  tllnl tnii'x , A 1 mtlx,  during 
nii'lt  Inn  mill  solidification. 


U«i  liquid  solid  Interface  Is  oqunl  to  (•';  • It  where  t.'g 
Is  I hi*  temperature  gradient  In  tho  liquid  id  tli«  liquid 
solid  Interface.  Tim  fuiu'tlon  (•' / • H i /'  (cooling  ride 
divided  by  square  of  absorlmd  heal  flux)  ex  fractional 
distance  solidified,  t /.  mux,  can  be  determined  from 
Kin.  It).  Initial  solidification  Interface  velocity  Is 
zero  but  increases  rapidly  with  fractional  distance 
solidified;  that  is,  on  an  expanded  « /.mux  scale,  the 
curves  in  Fig.  10  pass  through  the  origin.  As  noted 
from  the  temperature  profiles  presented  III  V'lgs.  3, 
ii.  and  7.  all  the  superheat  In  the  liquid  Is  lost  (t • / - It 


lie  comes  /.oro)  sometime  before  the  surface  layer  com 
plainly  solidifies.  Average  cooling  rate  In  the  liquid 
phase  at  any  time  during  solidification  rs  fractional 
distance  solidified  us  shown  In  Fig.  11.  This  cooling 
rate  Is  a maximum  at  the  beginning  of  solidification 
when  temperuture  gradients  In  the  liquid  are  high. 

Temperature  gradient  to  growth  rate  ratios  rs  (rue 
tlonul  melt  depth.  /,  /.mux  for  various  fractional  ills 
tames  solidified,  i /.,  are  shown  in  I'lg.  12  Kor  a 
given  melt  depth,  / /.  tm,x.  (>'/  It  is  a muximum  at  the 
beginning  of  solidification  and  drops  to  zero  sometime 


I'lg.  10—' The  ratio  of  cooling  rate  in  (lie  liquid  al  I lie  liquid 
solid  Interface  to  the  square  of  absorluxl  heal  liux,  (,'/  • N 
rs  fracllonal  distance  solidified,  « //  lmlx,  after  a maximum 
melt  depth  ts  achieved. 


pig,  1 1— The  ratio  of  average  cooling  rale  In  the  liquid  to  the 
square  ot  absorluxl  heat  llux.  (Mj  •ffWlrg/y*  i s fractional 
distance  solidified,  i / max- 
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Fig.  12— Temperature  gradient  to  growth  rate  ratio  at  the 
liquid-solid  interface,  G ^//< , as  fractional  molt  depth, 
l./L  max,  for  various  fractional  distances  solidified,  c/A. 

near  the  end  of  solidification  as  superheat  in  the 
liquid  is  lost.  Data  in  Fig.  12  also  show  that  O^/R  in- 
creases as  melt  depth  approaches  the  maximum  melt 
depth,  L /-max  1.0. 

Finally,  cooling  rates  for  solidification  of  noncrys- 
talline structures  were  calculated  by  setting  the  heat 
of  fusion,  II.  equal  to  zero  (see  Fig.  13).  Total  aver- 
age cooling  rate  was  obtained  by  averaging  over  melt 
depth  and  temperature  until  the  surface  reached  one- 
half  of  the  melting  temperature. 

VI.  DISCUSSION 

The  important  influence  of  high  cooling  rates  on  re- 
finement of  solidification  microstructure,  increased 
solid  solubility,  formation  of  nonequilibrium  crystalline 
and  noncrystalline  phases  has  now  been  well  estab- 
lished. Yet.  direct  correlation  of  temperature  profiles 
during  solidification  with  microstructure  has  been 
hartpered  due  to  experimental  difficulties.  The  availa- 
bility and  continuous  development  of  high  intensity 
power  sources  such  as  CO>  laser  may  alleviate  these 
problems  in  due  course.  First,  surface  melting  ex- 
periments carried  out  under  controlled  conditions 
would  lead  to  predictable  melt  depths  in  intimate  con- 
tact with  the  substrate  below.  These  are  two  important 
variables  that  cannot  be  controlled  in  splat  cooling. 
Second,  variation  in  heat  flux  would  permit  directional 
solidification  under  a variety  of  temperature  gradient 


and  growth  rate  conditions.  Finally,  isolation  of 
growth  from  nucleation  coupled  to  experimental  tech- 
niques yet  to  be  developed  for  accurate  temperature 
measurements  would  greatly  enhance  our  fundamental 
understanding  of  rapid  solidification  phenomena. 

The  results  presented  above  show  the  significant 
range  of  melting  and  solidification  variables  that  are 
available  with  variation  of  absorbed  heat  flux  when 
heat  flow  is  construed  to  be  unidirectional.  The  general 
trends  established  in  this  study  are  summarized  in 
Table  II. 

When  the  Absorbed  Heat  Flux  is  Increased  by 
an  Order  of  Magnitude,  the  Following  Effects  are 
Noted  on  the  Other  Variables: 

a)  Time.  The  various  times.  /,,,  and  fmax  de- 
crease by  two  orders  of  magnitude.  In  general,  the 
surface  of  a material  reaches  a given  temperature 
when  r/vT is  kept  constant.  The  significance  of  this 
finding  is  that  at  higher  power  inputs,  less  time  is 
available  for  diffusion  of  heat  into  the  metal  sub- 
strate-absorbed heat  is  concentrated  near  the  surface 
of  the  material  leading  to  steeper  temperature  gradi- 
ents. 


0 0.2  0.4  0.6  0.8  1.0 


Fig.  t:)— The  ratio  of  total  average  cooling  rate  to  the  square 
of  absorbed  heat  flux  rs  fractional  distance  melted,  A /A  max 
during  solidification  of  a noncrystalline  solid. 


Table  II  The  Effect  of  Change  in  Absorbed  Heat  Flux  on  Other  Variables 


Order  id  magnitude  increase  in  q results  in  the  following  changes 
if  rime  Meli  depth  H <> f (>/  • R (1//R 
lot  10*4  im  lot  lot  io’t 

f Increase 
1 deciease 
'*  ho  change 
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b>  MeH  Depth.  The  melt  depth  decreases  by  one 
order  of  magnitude  it  is  inversely  proportional  to  the 
absorbed  heat  flux.  As  an  example,  note  that  the  slope 
of  plots  of  maximum  melt  depth  vs  absorbed  heat  flux 

th'is  f °H  pU>t  *FiK'  11  18  1 The  Ke,ieral  forms  of 

this  finding  are  presented  in  Figs.  4 and  5.  Figure  4 

shows  that  for  a given  substrate  material,  the  di- 
mensionless melt  depth.  Hi.  „ during  melting 
and  subsequent  solidification  is  a constant  at  a given 
dimensionless  time.  I /,„.  Or.  for  a given  „v7.  the 
dimensionless  melt  depth.  /.  Z.max.  i8  a constant,  see 
r ig.  5. 

c)  Temperature  Gradient.  The  temperature  gradi- 
ents in  the  liquid  and  the  solid  regions  increase  bv  an 
order  of  magnitude.  This  fact  is  illustrated  in  Figs.  6 
7.  and  8 where  the  temperature  distributions  are  pre- 

nodal 'jxj s i Uon  °f  din,ension,e88  temperature  a.  each 

dl  lnterface  Velocity.  The  interface  velocity  both 
uring  melting  and  solidification  increase  by  an  order 
of  magnitude.  This  fact  can  be  directly  deduced  from 
the  above  and  is  illustrated  in  Fig.  9. 

e)  Cooling  Rates.  Cooling  rates  in  the  liquid  at  the 
lquid-solid  interface.  (.1,  ■ l{.  during  solidification 
increase  by  two  orders  of  magnitude,  see  Fig.  10. 
Average  cooling  rate  in  the  liquid  region  and  total 
average  cooling  rate  during  formation  of  crystalline 
and  noncrystalline  solids,  respectively,  follow  similar 
rends,  see  Figs.  11  and  13.  The  data  in  Fig.  13  show 
that  cooling  rates  at  a given  power  input  tend  toward 
a maximum  as  melt  depth  approaches  zero. 

f»  «.  The  ratio  of  temperature  gradient  in  the 
liquid  at  the  liquid-solid  interface  to  solidification  in- 
terface velocity  remains  constant. 

Some  experimental  evidence  f„  date  is  in  qualitative 
agreement  with  the  general  trends  noted  above  For 
example,  recent  work  by  Beck  vi  „/“  on  superalloy 
melts  produced  by  laser  irradiation  has  shown  that 
while  structures  near  the  root  of  a molten  zone  grow 
perpendicular  to  the  substrate  material,  solidification 
also  involves  growth  of  parallel  dendrites  near  the  ton 
of  this  zone.  This  observation  is  in  line  with  the 
predicted  loss  of  perpendicular  temperature  gradients 
in  the  liquid  sometime  during  solidification. 

The  effects  of  convection  in  the  liquid  and  heat  loss 
trom  the  melt  surface  were  also  investigated.  In- 
creasing the  liquid  metal  conductivity  by  a factor  of 
five  resulted  in  a similar  increase  in  maximum  melt 
depth,  whereas  incorporation  of  a heat  transfer  co- 
efficient at  the  liquid  air  interface,  an  order  of  mag- 
mtude  larger  than  that  calculated  for  radiation  had 
negligible  influence  on  the  temperature  gradients 
Finally,  kinetic  limitations  to  melting  and  solidifica- 
Bon  are  not  expected  to  significantly  effect  the  calcu- 
Uted  temperature  distributions,  if  we  assume  that 
both  melting  and  solidification  occur  in  the  con- 
tinuous regime  as  discussed  by  Cahn  «•/  ,,/17  then  the 
following  equation  is  applicable. 

,,  „ JDrHAT 

I’l 

where: 

A/  superheat  or  undercooling  at  the  interface 

during  melting  and  solidification,  respectively. 
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* = Holtz  man's  constant, 

No  3 Avagadro’s  number,  * 

<»  = step  height, 

d = parameter  that  relates  step  height  to  mean 
lump  distance. 

Maximum  undercooling  at  the  interface  would  thus 
e expected  at  the  highest  interface  velocities.  As- 
suming a large  absorbed  heat  flux  of  -10"’  W nr  a 
maximum  interface  velocity  of  -0.9  m s is  predicted 
from  Fig.  9 during  solidification  of  nickel.  The  cal- 
culated  undercooling,  for  ,1  - 00.  is  of  the  order  of 
' , ' *r°r  *arlJer  values  of  ,i  or  slower  interface 
velocities  the  corresponding  change  in  the  melting 
point  would  be  smaller.  h 


Vll.  SUMMARY 

A controlled  step-function  heat  input  via  a high  in- 
tensity continuous  CO,  laser  can  be  an  effective 
means  of  studying  rapid  solidification  of  thin  surface 
layers.  One  dimensional  heat  flow  calculations  show 
that  the  maximum  melt  depth  achieved  is  generally 
““  of  magnitude  less  than  the  heat  affected 
zone.  High  absorbed  heat  fluxes  lead  to  high  tempera- 
ture gradients,  shallow  melt  depths,  increased  inter- 

rre,0t;^f'  ?nd  si«niflcan»ly  higher  cooling  rates 
uring  solidification.  For  example,  the  surface  of  an 
aluminum  substrate  absorbing  a heat  flux  of  5 x 10" 

r®a.CheS  '!S  va»Joriza‘ion  temperature  in  approxi- 

denth ^Ji'rnV0  , S'  The  corresP°nding  maximum  melt 
pth.  solidification  interface  velocity,  o',  H and  maxi- 
mum average  cooling  rate  in  the  liquid  are  -65  pm  ~ 
0.45  nv  s.  3 x io"  K s m.  and  6.5  x 10"  K s.  respec- 
tively. For  the  same  power  input,  a melt  depth  of  ap- 

xTn-',111*1  14  2 is  achteved  in  approximately  2.42 

, 10  s;  the  corresponding  total  average  cooling  rate 
for  noncrystalline  solidification  of  this  melt  is  approxi- 
mately 7.5  x lo"  K s.  ppioxi 
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Ct>  specific  heat, .)  kg'1  K'1, 

distance  from  the  surface  to  the  last  nodal 
point  in  the  solid,  m or  tan, 
temperature  gradient.  K m'1, 

W.  temperature  gradient  in  the  liquid  at  the 

liquid-solid  interface,  K m'1, 

11  heat  of  fusion,  .1  kg'1, 

thermal  conductivity,  .1  s'1  m'1  K"1. 
melt  depth,  m or  pm, 

L'"a*  maximum  melt  depth  achieved  after  heat 
flux  is  removed  when  the  surface  reaches 
the  vaporization  temperature  of  the  mate- 
rial, m or  pm, 

L I-  max  dimensionless  melt  depth. 

PsNL  'll,,,  dimensionless  melt  depth. 

" increments  of  space,  dimensionless 

total  number  of  space  increments,  dimen- 
sionless, 

<1  absorbed  heat  flux,  w m*3, 

number  of  space  increments  in  the  liquid 
region,  dimensionless, 

W interface  velocity,  m s'1, 

distance  from  the  liquid-solid  interface  to 
the  last  nodal  point  in  the  solid,  m or  pm. 
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I ' total  time  after  application  of  heat  flux,  s, 

I/,,  time  for  the  surface  of  a material  to  reach 

its  melting  temperature,  s. 

/max  time  for  a material  to  achieve  a maximum 

melt  depth,  s. 

/,,  time  for  surface  of  a material  to  reach  its 

vaporization  temperature,  s. 

T temperature.  K. 

T„  ambient  temperature,  K, 

7\\/  melting  temperature,  K. 

Ts  surface  temperature,  K. 

T,,  vaporization  temperature.  K. 

.V  distance  from  the  surface,  m or  pm, 

i»  thermal  diffusivity  k/ p cy,  nr  s'1, 

t position  of  solidification  interface  meas- 

ured from  melt  depth.  /. , m or  pm, 
t*  I'  l'\i  temperature  excess  above  the  melting 
temperature,  K, 

fin  ni  temperature  excess  at  space  point  >1.  time 

point  m.  K. 

p density,  kg  m~\ 


Subscripts 

/ liquid  material, 

s solid  material. 
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